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SIMULTANEOUS EQUATION SOLUTION 



Author: D. N. Leeson 

WorLd Trade Europe Corporation 
3-5, Cite^du Re tiro 
Paris 8 e , France 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for the IBM Data Processing Systems. If such announce- 
ment indicates a change to the program decks or tapes, a 
complete new program, if needed, should be requested 
from the Program Distribution Center. 



DECK KEY 



L. Program Deck 

2. Sample Problem Deck 
(5 Cards) 



SIMULTANEOUS EQUATION SOLUTION (Card) 



Author: D. N. Lee son 

WorLd Trade Europe Corporation 
3-5, Cite" du Retiro 
Paris 8 e , France 

Direct Inquiries to: Author 

A. Purpose /De scr ipt ion : This program generates the soLutions to 
a linear system of maximum r ank 39 x 39. One may have 99 
constant vectors per matrix of coefficients. 

B. Method : Calculation of the product matrix. Arithmetic; floating. 

C. Restrictions and Range: 39 x 39 

Accuracy : Rounding error for very large systems noticeable. 

E. Machine Configuration : 1620 with 1622 attachment. Division 
feature not required. 

F. Program Requirements ; AH of core is required for the maximum 
problem. 

G. Source Language : This program uses SPS Language, and is 
non-relocatable. 

H. Program Execution Time: Variable dependent upon problem size. 

I. Check-Out Status: N/A 

J. Sample Problem Running Time: N/A 

K. Comments : This program and its documentation were written » 
by an IBM employee. It was developed for a specific purpose and 
submitted for general distribution to interested parties in the 
hope that it might prove helpful to other members of the data 
processing community. The program and its documentation are 
essentially in the author's original form. IBM serves as the 
distribution agency in supplying this program. Questions 
concerning the use of the program should be directed to the 
author ' s attention . 



DESCRIPTION 



This program generates solutions to a linear system of maximum 
matrix 39 X 39 with from 1 to 99 constant vectors. Division feature and 
indirect addressing are not necessary. The solution is accomplished 
without matrix inversion by the product matrix method. This method 
requires only m 3 multiplications and divisions as opposed to inversion 
which requires m 3 multiplications and divisions. The typed output of 
this program is the constant vectors and associated solution vectors. 
All input and output is in excess 50 floating point notation. 

CARD PREPARATION 

Data input consists of 1 header card,l or more matrix cards and 
I or more constant vector cards per constant vector set. 

HEADER CARD 

Cols. 1-2: xx - number if rows S number 

of columns in the matrix of coefficients. 

Cols. 3-4: xx - number of constant vectors 

MATRDC CARDS: The matrix of coefficients by rows in floating point 
form with high order digits flagged. Utilise as 
many columns and as many cards as needed with 
eight elements per card . A matrix of 6X6 will use 
4 full cards and 40 columns of a 5th. A matrix of 
2X2 uses only 1/2 of one card. A maximum matrix 
uses 190 full cards and 10 columns of the I9lst. (39 X 39) 

-I- 



CONSTANT VECTOR CARD(S): The constant vector in floating point form 

with high order digits flagged. Utilize as 
many columns and as many cards as needed 
with eight elements per card. A matrix of 
16 X 16 with 1 constant vector of 6 elements 
requires 2 cards for constant vector 
description. 

If more than 1 constant vector exists 
begin a new set of constant vector card{s) 
for each constant vector 



EX 



2X J. 3Y ~ 
9X - 2Y 



«0 



CARD 1 0203 

2 5120000000513000000051900000005120000000 

3 51700000005140000000 

4 51600000005212000000 

5 51300000005180000000 



-2- 



OPERATION: 



1. Clear memory 

2. Load program deck 

3. Program halts after loading. 

4. Ready data in card- read hopper 

5. Press start. 

6. Program solves for all constant vectors and halts. If 
further systems are to be solved, ready reader with 
data and press start. 
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